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Abstract: We report the results of a successful 7 hour 1.4 GHz VLBI experiment using two new 
stations, ASKAP-29 located in Western Australia and WARK12M located on the North Island of 
New Zealand. This was the first geodetic VLBI observing session with the participation of these new 
stations. We have determined the positions of ASKAP-29 and WARK12M. Random errors on position 
estimates are 150-200 mm for the vertical component and 40-50 mm for the horizontal component. 
Systematic errors caused by the unmodeled ionosphere path delay may reach 1.3 m for the vertical 
component. 
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1 Introduction 



As part of the Australian and New Zealand joint 
bid to host the multi-billion dollar Square Kilo- 
metre Array (SKA), both coimtries are investing 
heavily in advanced technologies for radio astron- 
omy. In Australia this follows a strong tradi- 
tion in radio astronomy and is expressed in the 
construction of the Australian SKA Pathfinder 
(ASKAP) on the western edge of the Australian 
continent (Johnston et al. 2008). Radio astron- 
omy in New Zealand has links stretching back to 
the work of Elizabeth Alexander on solar radio 
emission (Alexander 1946). John Bolton and Gor- 
don Stanley used a cliff interferometer to obtain 
rising and setting records of various radio sources 
(Bolton 1982); measurements made near Sydney 
and Auckland allowed them to identify mysterious 
"radiostars" with well known supernova remnants 
and galaxies (Bolton, Stanley, & Slee 1949). The 
first successful VLBI experiment between Aus- 
tralia and New Zealand was made in 2005 with a 
6-m radio telescope near Auckland (Karaka) and 
the Australia Telescope Compact Array (ATCA) 
(Gulyaev et al. 2005; Tingay et al. 2006). In 2008, 
Auckland University of Technology commissioned 
a new 12-m antenna at Warkworth near Auckland, 
— the New Zealand's first research capable radio 



telescope (Gulyaev & Natusch 2009). 

A joint Australia and New Zealand collabo- 
rative project is focused on developing both the 
ASKAP and Warkworth facilities as part of re- 
gional and global Very Long Baseline Interferom- 
etry (VLBI) arrays, for the purposes of both as- 
tronomy and geodesy. Due to the dearth of land 
mass in the Southern Hemisphere, both astronom- 
ical and geodetic VLBI has suffered in the past, 
with the radio telescopes used for VLBI restricted 
to the south east "corner" of the Australian conti- 
nent, one telescope in South Africa and occasional 
use of radio telescopes in Antarctica and South 
America. 

The additional capability gained by adding ra- 
dio telescopes located in Western Australia and 
New Zealand is substantial, increasing the angu- 
lar resolution of the Australasian array by a factor 
of approximately four, admitting a range of as- 
tronomy science goals described in Johnston et al. 
(2008). 

High precision astrometry is likely to be an im- 
portant part of the science case for the high angu- 
lar resolution component of the SKA, in particular 
for precise determination of the distances of radio 
pulsars in the galaxy, to be used in various tests of 
fundamental physics (Smits et al, 2010). In order 
to achieve the astrometric performance required. 
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the SKA will need to undertake astrometry util- 
ising distributed elusters of small antennas oper- 
ating as phased arrays. This is a departure from 
standard precision astrometry and our work using 
ASKAP as part of a VLBI array for astrometry 
will be an important testbed for the demonstra- 
tion of these techniques. 

The 12-m radio telescope WARK12M^ is in- 
tended to be used as a part of the Australian Long 
Baseline Array (LEA), for spacecraft monitoring, 
and for VLBI observations in the framework of 
the International VLBI Service for Geodesy & As- 
trometry (IVS)^ and the AuScope project'^. 

ASKAP-29 is an element of the ASKAP array of 
36 identical dishes using Phased Array Feed (PAF) 
technologies, expected to be fully operational by 
2013. ASKAP will undertake very wide field sur- 
vey science in continuum and spectral line modes 
and is also intended to be used as a part of re- 
gional and global VLBI networks for a variety of 
projects''. 

As a part of commissioning the new antennas, 
the positions of antenna reference points should 
be determined. A reference point is defined as the 
point of the projection of the movable elevation 
axis onto the fixed azimutlial axis. For the analy- 
sis of VLBI source imaging experiments made in a 
phase-referencing mode, the projection of a base- 
line vector, (i.e. vector between antenna reference 
points) to the tangential image plane should be 
known with errors not exceeding tens of centime- 
tres, otherwise the image will be smeared Chariot 
et al. (2002). For astrometry applications, angu- 
lar position accuracies of tens of fias arc required, 
making the requirements on the accuracy of sta- 
tion positions much more stringent: 5-10 mm. 
One way of estimating the position of the antenna 
reference points is through analysis of a combi- 
nation of a ground survey of markers attached to 
the antennas from a local network around the sta- 
tion and GPS observations from the points at the 
local network (see e.g., Sarti, Sillard & Vittuari 
2004; Sarti, Abbondanza & Vittuari 2009). An- 
other way to estimate station positions is to use 
the VLBI technique itself to determine group de- 
lays and then derive reference point positions from 
these group delay measurements. The advantage 
of this approach is that it also provides useful di- 
agnostics on the VLBI equipment. 

The first fringe test experiment between ASKAP- 
29 and WARK12M was made in April 2010. First 



-'^wark12m and ASKAP-29 are identifiers for specific 
VLBI antennas: near Warkworth, New Zealand and at the 
Murchison Radio- Astronomy Observatory, Western Aus- 
tralia, respectively. 

^http : //ivscc . gsf c .nasa. gov/ 

^http : //www . auscope . org. au/ 

''http : //www. atnf . csiro . au/SKA/ 
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Table 1 Specifications of an ASKAP antenna 



Antenna Type 

Reflector 
Surface accuracy 
Mount 

Height of elev axis 

Azimuth axis range 
Azimuth axis speed 
Azimuth axis accel. 

Elevation axis range 
Elevation axis speed 
Elevation axis accel. 

Polarization axis range 
Polarization axis speed 
Polarization axis accel. 



fuUy-steerable, 
prime focus 
12-metre paraboloid 
0.5 mm rms or better 
3 axes: polarization, 
elevation, azimuth 
7.51 m above 
concrete foundation 

-180° to +360° 

3°/s 

3%^ 

15° to 89° 
l°/s 

1°/S2 

-180° to -H80° 
3% 

3°/s^ 



fringes on baselines between ASKAP-29 and MOPRA 
were found on April 22, 2010 and on baselines 
to WARK12M on the following day. This success 
prompted three first science experiments using the 
full Long Baseline Array (LBA) network with the 
participation of the two new stations: 1) imaging 
observations of PKS 1934-638 on 29-30th April 
(Tzioumis et al. 2010a), 2) a geodetic experiment 
on May 07, 2010, and 3) imaging onbservations of 
Cen-A on May 09, 2010 (Tingay et al, 2010, in 
preparation). 

We report here results from the geodetic exper- 
iment that was conducted at ASKAP-29 and dur- 
ing the first geodetic observing session in May 07, 
2010. The goal of this experiment was to deter- 
mine the position of the antennas with decimeter 
accuracy and to collect diagnostic data. The char- 
acteristics of new antennas arc presented in section 
2. The experiment and its analysis are described 
in sections 3 and 4. Concluding remarks are made 
in section 5. 

2 New Antennas 

2.1 ASKAP-29 

The Austrahan SKA Pathfinder (ASKAP) is an 
advanced technology radio telescope currently un- 
der construction at the Murchison Radio Observa- 
tory (MRO) in Western Austraha. The ASKAP 
array comprises 36 fuUy-steerable dish antennas of 
12 metre aperture fitted with phased array feed 
(PAF) receivers, giving an instantaneous field of 
view of approximately 30 square degrees with an 
instantaneous bandwidth of 300 MHz in the range 
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700-1800 MHz. 

The first antenna, ASKAP-29, was commis- 
sioned in March 2010 and was temporarily fitted 
with a conventional horn (single-pixel) L band feed 
(1.4 GHz) and uncooled receiver. Five additional 
antennas were completed in Q4 2010 with the bal- 
ance to be commissioned throughout 2011. 

The antenna is of a prime focus design with an 
unshaped paraboloidal reflector of f/D 0.5. The 
original specifications on surface accuracy were 
1 mm root mean square (rms). However, mea- 
surements showed that the surface accuracy of the 
ASKAP-29 reflector and six other identical antennas 
manufactured by November 2010 is better than 
0.5 mm rms over the full elevation range. There is 
some flexure in the position of the focus at a level 
of few mm over the full elevation range, but this 
can be compensated for if required. This level of 
surface accuracy implies that the antenna will be 
usable at frequencies as high as Ka band (30 GHz) 
and possibly at Q-band (43 GHz). A quadripod 
supports the prime focus platform with maximum 
payload of about 200 kg. Additional bracing sup- 
ports on each leg have been added to the original 
design to reduce flexure. 

An unusual feature of the antenna is the three- 
axis mount. Above the familiar azimuth and el- 
evation axes is the so-called polarization axis, al- 
lowing the entire reflector to be rotated about its 
optical axis. This design was specified to allow full 
parallactic angle tracking for high dynamic range 
in imaging observations, as it improves modeling 
of far-out sidelobes caused by the quadripod. An 
equatorial ( "HA-Dec" ) design naturally maintains 
a constant orientation of its focal plane with re- 
spect to the celestial sphere when tracking a fixed 
point in celestial coordinates. However this type of 
mount usually suffers restricted sky coverage due 
to limited rotation in hour angle. 

Another feature of the antenna is the relatively 
large pedestal specified to accommodate the ana- 
logue and digital electronics for the PAF receiver 
with its 188 individual ports. The 54th Research 
Institute of China Electronics Technology Group 
Corporation (known as CETC54) are contracted 
to design, manufacture and install all thirty-six 
antennas at the MRO. 

The baseband data at ASKAP-29 were recorded 
directly to disk utilising a custom recording sys- 
tem developed at Curtin University. The record- 
ing systems were assembled from consumer off- 
the-shelf components. The recorders were inter- 
faced with the CSIRO supplied analogue down- 
converters. The systems recorded the data to 
disk at the native sampler/digitizer sample size of 
16 bits. 

The 16 bit sample size gives an aggregate data 
rate when sampled at the Nyquist-Shannon rate 




Figure 1 ASKAP antenna. 

over 2 polarization channels of 4 Gbit/s. The 
recorder system comprised a Signatec PX14400 
PCIc sampler/digitizer mounted in a server class 
computer. The systems were built in a 16 disk 
chassis populated with SATA hard drives. The 
recording system utilised a Linux software RAID, 
and the 16 disks where subdivided into 8 disk raid 
units. The recorded data were converted to the 
LBA format with 2 bit samples after the exper- 
iment using software employing automatic gain 
control in the conversion process in order to match 
the configuration of LBADR recorders used at 
other stations. 

2.2 WARK12M 

The New Zealand 12-m radio telescope (see Fig- 
ure 2) is located some 60 km north of the city of 
Auckland, near the township of Warkworth. It was 
manufactured by Patriot Antenna Systems (now 
Cobham Antenna Systems), USA. The antenna 
specifications are provided in Table 2. The ra- 
dio telescope was originally designed to operate at 
S (2.3 GHz) and X bands (8.4 GHz) and it was 
supplied with an S/X dual-band dual-polarisation 
feed. Photogrammetry observations showed that 
the surface rms accuracy is 0.35 mm. Therefore, 
the antenna is suitable for observations at frequen- 
cies as high as 43 GHz. It is equipped with a digi- 
tal base band converter (DBBC) developed by the 
Italian Institute of Radio Astronomy, Symmetri- 
com Active Hydrogen Maser MHM-2010 (75001- 
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Table 2 Specifications of an wark12m antenna 



Antenna type 


Fully-stecrable, 




dual-shaped Cassegrain 


Manufacturer 


Cobham/Patriot, USA 


Main dish Diam. 


12.1 m 


Secondary refl. Diam. 


1.8 m 


Focal length 


4.538 m 


Surface accuracy 


0.35 mm 


Pointing accuracy 


18" 


Frequency range 


1.4—43 GHz 


Mount 


alt- azimuth 


Azimuth axis range 


90° ± 270° 


Elevation axis range 


4.5° to 88° 


Azimuth axis max speed 


5°/s 


Elevation axis max speed 


l°/s 


Main dish F/D ratio: 


0.375 



114) and with the Mark5B-|- data recorder de- 
veloped at MIT Haystack Observatory. However 
these were not installed at the time of this experi- 
ment, so data were recorded with the same equip- 
ment as installed at ASKAP-29. 

The antenna elevation axis is at a height of ap- 
proximately 7.1 m above the ground level. The 
elevation axis is supported by a pedestal that is 
of steel construction. The pedestal is essentially 
a steel cylinder of ~ 2.5 m diameter. Apart from 
the pedestal all other components of the antenna 
(the reflector and feed support structure) arc con- 
structed of aluminum. The support foundation 
for the antenna is a reinforced concrete pad that 
is 1.22 m thick by 6.7 x 6.7 meters square. The 
ground that the foundation is laid on consists of 
weathered sandstone/mudstone, i.e. is of sedimen- 
tary origin, laid down in the Miocene period some 
20 million years ago. 

The radio telescope is directly connected to 
the regional advanced network KAREN (Kiwi Ad- 
vanced Research and Education Network), which 
provides a fast connectivity between New Zealand's 
educational and research institutions (Weston, Na- 
tusch & Gulyaev 2010). At the time of the experi- 
ment, the Trans- Tasman (New Zealand- Australia) 
connectivity from the observatory was 155 Mbps. 
By the end of 2010 this was upgraded to 1 Gbps. 

A preliminary survey has been conducted in 
collaboration with the New Zealand Crown Re- 
search Institute (CRI), GNS Science and Land In- 
formation New Zealand (LINZ) to determine an 
initial estimate of the reference point of the VLBI 
site WARK12M. This reference point is defined as 
the intersection of the azimuth and elevation axes 
of the telescope. A real-time kinematic (RTK) 
GPS method was used to derive the position with 
respect to the GPS station WARK. 

WARK was established in November 2008 at 



Figure 2 wark12m antenna. 



the radio telescope site and is one of thirty nine 
PositioNZ network stations (Bhck & CoUett 2010) 
in New Zealand^. All data received from the Posi- 
tioNZ stations are compiled into 24 hour sessions 
and are processed to produce daily positions for 
each station in terms of ITRF2000^. The coordi- 
nates for WARK used in the following calculations 
were derived by averaging the daily coordinate so- 
lutions for February 19 through March 09, 2010. 

The RTK reference receiver was set up in an 
arbitrary location with clear sky view and was con- 
figured to record raw observations in addition to 
transmitting real-time corrections. This station 
was later post-processed with respect to WARK 
and all RTK rover surveyed positions were subse- 
quently adjusted relative to the updated reference 
position. 

Several points on the rim of the main reflector 
were identified and each point was measured sev- 
eral times with the RTK rover while the telescope 
was repositioned in elevation and azimuth between 
each measurement. The rover GPS antenna was 
mounted on a 0.5 m survey pole and was hand 



^http: / /www. llnz. govt. nz/geodctic/positionz/ 
"^For more details sec 

http : //geonet . org .nz/resources/gps/ 

gps-processing-notes .html 
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held for each measurement. Access to the rim of 
the reflector was achieved with a hydrauHc cherry 
picker. 

The sequence of observation for determination 
of the horizontal axis was as follows: The tele- 
scope azimuth axis was held fixed (nominally 0°). 
A point near the highest edge of the reflector was 
identified and measured with the telescope in 4 po- 
sitions of elevation, from almost zenith (^88°) to 
as high as the cherry picker could reach 38°). 
This was repeated a second time with a point iden- 
tified on the edge of the reflector to one side of the 
telescope. Five positions of elevation were mea- 
sured at this point, from 10° to 80°. 

The sequence of observation for determina- 
tion of the vertical axis was as follows: The tele- 
scope elevation axis was held fixed (nominally 
80°). Three points around the edge of the refiector 
were identified. The telescope azimuth axis was 
rotated into three positions such that each iden- 
tified point could be measured consecutively from 
the surveyor's location in the cherry picker. The 
cherry picker was then repositioned twice around 
the perimeter of the telescope and the measure- 
ments were repeated at each cherry picker loca- 
tion. This provided three measurements of each 
identified point with a fixed telescope zenith and 
varying azimuths. 

The resulting points from these measurements 
describe two vertical circles of rotation which de- 
fine the movable elevation axis and three horizon- 
tal circles of rotation which define the fixed az- 
imuthal axis. The coordinates for all subsequent 
calculations were retained as geocentric Cartesian 
coordinates to avoid any possibility of errors re- 
lated to transformation of projection. 

The following method was used to determine 
the axes and their intersection point from these ob- 
servations. The equation of a circle from 3 points 
(P. Bourke, 1990, internal memo^) was used to cal- 
culate all possible combinations of three observed 
points which define a circle of rotation. A simple 
mean was taken for all horizontal axis definitions 
and all vertical axis definitions. The mid point of 
the closest point of approach of each axis with the 
other was used as the final estimate of the point of 
intersection. The distance between closest point of 
approach on each axis was calculated to be 24 mm. 
Based on the variation of results for different com- 
binations of survey points, we estimate that the 
accuracy of the determined intersection point is 
within 0.1 m. 

In summary, the following coordinates of the 
intersection of the azimuth and elevation axes for 
the radio telescope WARK12M were derived in terms 
of ITRF2000 at the epoch of the survey (March 

^Available at http : //local .wasp. uwa. edu. au/'pbourke/ 
geometry/circlef rom3 



2010): 

X = -5115324.5 ± 0.1m 

Y = 477843.3 ± 0.1m 

Z = -3767193.0 ± 0.1m 

It is intended that the radio telescope refer- 
ence point coordinates will subsequently be re- 
determined to a higher accuracy with the use of a 
variety of terrestrial and GNSS survey techniques 
Dawson, Johnston & Twilley (2004) and a more 
rigorous least squares analysis of the observations. 
Four geodetic survey monuments have been built 
within 15-20 m from the antenna pedestal for this 
purpose. 



3 Geodetic VLBI Experiment 

Seven stations participated in the 7 hour long 
VLBI experiment vtl4a on May 07, 2010, with 
the goal of determining station coordinates. The 
stations are listed in Table 3 and shown in Fig- 
ure 3. Usually geodetic experiments are made in 
two bands recorded simultaneously at 2.2-2.3 GHz 
and 8.0-8.8 GHz. Since group delay errors are 
reciprocal to the frequency range, a wide range, 
720 MHz, is used in VLBI experiments dedicated 
to geodesy. Radio waves are delayed in the iono- 
sphere, and the magnitude of this delay depends on 
total electron content along the paths — a highly 
variable quantity. Since the contribution of the 
ionosphere to path delay is reciprocal to the square 
of effective frequency, simultaneous observations 
at two widely separated frequencies allow the for- 
mation of ionosphere-free combinations of group 
delays. 

However for this experiment, the data were 
recorded in the band 1368-1432 MHz in dual po- 
larisation (i.e. a single 64 MHz channel) . Observa- 
tions within a narrow-band at a low frequency are 
highly unusual for geodetic style experiments since 
they are not able to provide accuracy comparable 
with dual-frequency, wide-frequency range experi- 
ments. However, we ran this experiment because 
in May 2010 the new antennas had only L band re- 
ceivers, and this experiment provided an opportu- 
nity not only to make a coarse estimate of station 
coordinates, but also to test the equipment. The 
single frequency setup is a limitation of the VLBI 
baekend available at both ASKAP-29 and WARK12M 
at the time of observations. 

The experiment was scheduled using sur_sked 
software program in the VLBI geodetic mode. A 
pool of 113 sources with correlated flux densities 
at X band greater than 700 mJy based on results 
of the ongoing LBA Calibrator Survey (Petrov 
et al. 2011a) observing campaign was used. The 
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Table 3 The network stations 



Name 


(f>gd 


A 


Diam 


ASKAP-29 


-26°. 69 


116°.64 


12 


m 


ATCAPN5 


-30°. 31 


149°.56 


22 


m 


HOBART26 


-42°.81 


147°.44 


26 


m 


MOPRA 


-31°.27 


149°.10 


22 


m 


PARKES 


-33°.00 


148°.26 


64 


m 


WARK12M 


-36°.43 


174°.66 


12 


m 



scheduling goal was to have a uniform sky cover- 
age at each station. This was achieved by putting 
a next source in the schedule that has a maximum 
score among other candidate sources. The score 
was assigned according to 



S=- 



0.66 



(1) 



where ts is maximum slewing time for all partici- 
pating antennas in seconds, Dmin is the distance 
in radians to any source observed in the past. If a 
source was observed during the previous 150 min- 
utes, it was barred from consideration as a candi- 
date. Each scan had a nominal duration of 150 s. 
Every hour a set of 4 atmosphere calibration scans 
was observed: two scans with all antennas with 
their elevation in the range [12°, 45°], one scan 
at elevations [32°, 45°] and one scan at elevations 
[45°, 90°]. The scheduling algorithm for each set 
found all combinations of calibrator sources that 
fell in the elevation ranges and selected the se- 
quence of four objects that minimized the slewing 
time. The purpose of these calibration observa- 
tions was 1) to serve as amplitude and bandpass 
calibrators; 2) to improve robustness of estimates 
of the path delay in the neutral atmosphere. 



^20 




>0 



Figure 3 The LBA network used in vtl4a experiment. 



4 Data Analysis 

4.1 Data transfer and correlation 

Following the observations, the data from ATCAPN5, 
HOBART26, MOPRA and PARKES were transferred, 
with the support of ARCS (the Australian Re- 
search Collaboration Service), via fast network 
connections to a petabyte data store operated 
by iVEC in Western Australia. The data from 
WARK12M were transferred from the telescope us- 
ing the trans- Tasman link directly to Perth. The 
data from ASKAP-29 were transported to Perth by 
car, as the fast network to the MRO at the time of 
the experiment was under construction. All data 
were later transferred to disks mounted on the 
correlation cluster at Curtin University and cor- 
related using the DiFX software correlator (Deller 
et al. 2007). Initial antenna positions for WARK12M 
and ASKAP-29 determined from GPS receivers were 
used for correlation. Other station positions are 
known from VLBl solutions with accuracy better 
than one centimeter (Petrov et al. 2009a). The 
data were correlated using an integration period 
of 0.25 seconds and 512 frequency channels across 
the 64 MHz band. All four combinations of right 
and left circular polarization signals were corre- 
lated. 

4.2 Data quality control made at 
the correlator 

Initial data quality checking and clock search- 
ing was performed at the correlator using an au- 
tomated pipeline script implemented using the 
ParselTongue AIPS interface (Kettenis et al. 
2006). Hydrogen masers were used to provide 
the frequency standard at all stations except at 
ASKAP-29 where a lower stability (but still ac- 
ceptable) Rubidium frequency standard was used. 
This resulted in a noticeable correlation loss 
on baselines to ASKAP-29. All stations, except 
WARK12M observed with dual circular polariza- 
tions, while WARK12M observed dual linear polar- 
izations due to a malfunction of the linear to circu- 
lar conversion module, resulting in a reduced SNR 
on WARK12M baselines: ^ less sensitive than 
expected from a priori SEED estimates. ATCA 
participated as a phased array comprising 5x22-m 
dishes. Polarization isolation at parkes receiver 
did not work correctly during the first 2 hours of 
the experiment, but this was fixed after the linear 
to circular hybrid module was reset. 

4.3 Post-correlator data analysis 

The spectrum of cross correlation and autocorrela- 
tion for each scan and each baseline computed by 
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the correlator at a uniform two-dimensional grid 
of accumulation periods and frequencies was used 
for further processing. 

The fringe fitting procedure searches for phase 
delay Tp, phase delay rate fp, group delay Tg, and 
its time derivative fg that correct their a priori 
values used by the correlator model in such a way 
that the coherent sum of weighted complex cross- 
correlation samples over a given baseline and a 
given scan, c,;j 



Tn = 



2^ 2^ C^J X 



« 3 (2 

gl(woTp + "O-rp(ti-to) + (Wj-CJo)T"g + ( W j - ) Tg (* i - * ) 



reaches the maximum amplitude. Index i runs 
over time and index j runs over frequencies, wq 
and ^0 denote the angular reference frequency 
within the band and the reference time within a 
scan and Wij are weights that are defined as a frac- 
tion of processed samples in each accumulation pe- 
riod. Software VTM.A was used for amplitude cal- 
ibration and evaluation of phase and group delay, 
as well as their time derivatives. The algorithm 
implemented in VTM.A is described in detail in 
Petrov et al. (2011b). 

After computing phase and group delays and 
their derivatives, the spectrum of cross-correlations 
can be averaged over time or frequency. This pro- 
vides valuable diagnostics. The amplitude spec- 
trum at new stations ASKAP-29 and WARK12M 
shows a significant drop of the amplitude near the 
edges of the band as it is seen in Figure 4. The 
SNR losses due to the shape of the filter is at a level 
of 13%. A portion of the band with the fringe am- 
plitude less than 1 /4 of the maximum was filtered 
out in the final processing. 

Since the observations were made during low 
solar activity, the ionosphere did not cause sig- 
nificant decorrelation. Instability of the Rubid- 
ium frequency standard is clearly seen in Figure 5, 
and it caused a de-correlation on baselines with 
ASKAP-29 at a level of 4-5%. 

At the first step of geodetic data analysis, 
the spectrum of the cross-correlation function was 
screened. The edges of the bandpass and several 
spurious peaks near the center of the bandpass 
were filtered out. The complex bandpasses were 
evaluated and the fringe-fitting procedure was re- 
peated with the the refined bandpasses and with 
the bandpass mask applied. This procedure is 
described in full detail in Petrov et al. (2011b). 
Analysis of the amplitudes at all combinations 
of right and left circular polarizations, RR, LL, 
LR and RL, confirmed that the polarization setup 
at PARKES for the first two hours was incorrect. 
But it was also found that polarization impu- 
rity defined as ^\RLY + \LR\^I^\RRY + \LLY 
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Figure 4 The amplitude (upper plot) and phase (lower 
plot) of the cross correlation spectrum in a scan of source 
0537—441 at baseline ASKAP-29/parkes averaged over time 
after applying fringe search. Integration time 0.25 s 



was significant on baselines with other stations. 
For example, the polarization impurity at base- 
line ATCAPN5/MOPRA was at a level 7%, on base- 
lines ATCAPN5/HOBART2626, ATCAPN5/ASKAP-29 
and ATCAPN5/PARKES was 25-30%, and the base- 
line ATCAPN5/WARK12M was 100%. There- 
fore, the coherent combining dual-polarization im- 
proved the SNR at some baselines and degraded 
at others. For this reason, the datasets of group 
delays with RR and LL polarizations were used as 
independent experiments in further processing. 

The next step of the analysis pipeline is to com- 
pute theoretical path delays and form small dif- 
ferences between them and derived group delays. 
Computation of theoretical path delays was made 
with the use of the VTD software program^ that 
incorporates the state-of-the art geophysical mod- 
els. Our computation of theoretical time delays in 
general follows the approach of Sovers, Fanselow & 
Jacobs (1998) with some refinements. The most 
significant ones are the following. The advanced 
expression for time delay derived by Kopeikin & 
Schafer (1999) in the framework of general rel- 
ativity was used. The displacements caused by 
the Earth's tides were computed using the numer- 
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Table 4 Coordinates of ASKAP-29 and wark12m derived from analysis of experiment vtl4a. 



Station 




Station coordinates (m) 




X 


Y Z 


ASKAP-29 


-2556741.89 ±0.09 


5097441.23 ±0.14 -2847748.34 ± 0.11 


WARK12M 


-5115325.55 ±0.10 


477842.95 ± 0.05 -3767194.41 ± 0.09 



n 



-1 




_l I I l_ 



_l I l_ 
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Figure 5 Tlie phase of the cross correlation spectrum in 
a scan of source 0537—441 at baseline WARKl2M/pAaKES 
(upper plot) and ASKAP-29/parkes (lower plot) averaged 
over frequency after applying fringe search solutions. 



ical values of the generalized Love numbers pre- 
sented by Mathews (2001) following a rigorous al- 
gorithm described by Petrov & Ma (2003) with a 
truncation at a level of 0.05 mm. The displace- 
ments caused by ocean loading were computed 
by convolving the Greens' functions with ocean 
tide models. The GOT99.2 model of diurnal and 
semi-diurnal ocean tides (Ray 1999), the NA099 
model (Matsumoto, Takanezawa & Ooe 2000) of 
ocean zonal tides, the equilibrium model (Petrov 
& Ma 2003) of the pole tide, and the tide with 
period of 18.6 years were used. Station displace- 
ments caused by the atmospheric pressure loading 
were computed by convolving the Greens' func- 
tions that describe the elastic properties of the 
Earth (Farrell 1972) with the output of the atmo- 
sphere NCEP Reanalysis numerical model (Kalnay 



et al. 1996). The algorithm of computations is 
described in full detail in Petrov & Boy (2004). 
The displacements due to loading caused by varia- 
tions of soil moisture and snow cover in accordance 
with GLDAS Noah model (RodcU 2004) with a 
resolution 0.25° x 0.25° were computed using the 
same technique as the atmospheric pressure load- 
ing. The empirical model of harmonic variations 
in the Earth orientation parameters heo_20101111 
derived from VLBI observations according to the 
method proposed by Petrov (2007) was used. The 
time series of UTl and polar motion derived by the 
NASA Goddard Space Flight Center operational 
VLBI solutions were used a priori. 

The a priori path delays in the neutral atmo- 
sphere in the direction of observed sources were 
computed by numerical integration of differen- 
tial equations of wave propagation through the 
heterogeneous media. The four-dimensional field 
of the refractivity index distribution was com- 
puted using the atmospheric pressure, air temper- 
ature and specific humidity taken from the out- 
put of the Modern Era Retrospective- Analysis for 
Research and Applications (MERRA) (Schubert 
et al. 2008). That model presents the atmospheric 
parameters at a grid l/2°x2/3°x6'' at 72 pressure 
levels. 

In the initial least square (LSQ) solution, po- 
sitions of all stations, except parkes, were esti- 
mated, as well as coefficients of the expansion of 
clock function and residual atmosphere path delay 
in zenith direction into the B-spline basis of the 1st 
degree. The quality check revealed a clock break at 
station HOBART26. During the preliminary phase 
of the data analysis, outliers were eliminated and 
the baseline-dependent corrections to the a priori 
weights defined to be reciprocal to formal uncer- 
tainties of group delays were determined in such a 
way that the ratio of the weighted sum of squares 
of residuals to their mathematical expectation was 
close to unity. Group delays from RR and LL po- 
larization data were processed independently. 

The final LSQ solution used all VLBI group 
delays collected from 1984 through 2010^, in to- 
tal, 7.5 million values, including 824 group delays 
from this experiment. Positions of all stations, all 
sources, the Earth orientation parameters as well 



^Avaible at the IVS Data Center at 
http : // ivscc . gsfc .nasa. gov/products-dat a/index . html 
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as over 1 million nuisance parameters were esti- 
mated in a single LSQ run. Minimal constraints 
were imposed to require that net-translation and 
net-rotation over new position estimates of 48 sta- 
tions with long history with respect to positions 
of observations to these stations in the ITRF2000 
catalogue (Altamimi, Sillard & Boucher 2002) to 
be zero. This ensures that positions of all stations, 
including ASKAP-29 and WARK12M be consistent 
with the ITRF2000 catalogue. More details about 
parameter estimation technique can be found in 
Petrov et al. (2009b). 

The estimates of coordinates of ASKAP-29 and 
WARK12M on epoch 2010.05.07 are given in ta- 
ble 4. The errors reported in the table are the 
formal uncertainties from the LSQ solution com- 
puted in accordance with the error propagation 
law. Since there was only one experiment avail- 
able with rather an unusual setup, it is difficult to 
provide a realistic estimate of errors. 

The largest source of systematic error is the 
path delay in the ionosphere. Our attempt to use 
maps of the total electron content from GPS using 
the data product from the analysis center CODE 
did not improve the fit and shifted estimates of sta- 
tion coordinates at a fraction of the formal uncer- 
tainty (0.2cr). Analysis of dual-band VLBI exper- 
iments showed that global TEC maps above Aus- 
tralia and New Zealand provided by the Interna- 
tional Global Navigation Satellite System Service 
are not reliable, and the path delay computed from 
this model poorly represents the true delay in the 
ionosphere, at least during solar minimum (Petrov 
et al. 2011a). 

The adjustments to the residual zenith path 
delay in the atmosphere were in the range of 1- 
5 ns, while they arc typically in the range 0.03- 
0.1 ns for dual-band observations. Unlike the esti- 
mation of the mismodeled troposphere path delay, 
estimation of unmodeled ionosphere path delay in 
zenith direction from the observations themselves 
does not adequately represent the true path de- 
lay, since when antennas point at different direc- 
tions, the ionosphere piercing points with a typical 
height ~ 450 km may be located at distances up 
to 1000 km. 

For evaluation of the robustness of our coordi- 
nate estimates we made a trial solution and esti- 
mated the position of station ATCAPN5 indepen- 
dently. Estimates of the of ATCAPN5 position ap- 
peared within 2cr of the CATN5 pad. The position 
of CATN5 was derived from the position of the pad 
CATW104 determined from previous VLBI obser- 
vations and from the results of a local survey that 
measured coordinates of CATN5 with respect of 
CATW104 (Petrov et al. 2009a). This gave us a 
hint that the reported uncertainty may be underes- 
timated by a factor of 2, especially for the vertical 



component. In our final solution wc applied the ec- 
centricity vector between CATN5 and CATW104 
from the local survey. 

The data from RR and LL polarizations have 
independent random errors, but they share the 
same systematic errors. In another trial solution 
we processed RR and LL data separately. The po- 
sition differences of these two solutions were within 
la for horizontal component and 2a, or 0.4 m for 
the vertical component. 

However, the latter test is not sensitive to 
common systematic errors caused by unmodeled 
ionosphere path delays. In order to evaluate 
the magnitude of the ionosphere path delay er- 
rors, wc used prior observations at the baseline 
parkes/hobart26 under the IVS geodetic pro- 
gram. We selected 21 twenty four hour experi- 
ments in 2005-2010 — the period of solar min- 
imum. We ran three sets of 21 trial solutions. 
Within each set we used all experiments from 
1980 through 2010 with station HOBART26 re- 
moved from all but one experiment. The only 
experiment with hob art 26 was different in each 
solution within a set. In such an experiment all 
other stations, except HOBART26 and PARKES were 
removed. This setup emulated the case of deter- 
mining station coordinates from one experiment 
only. The reference set A used an ionosphere-free 
combination of X and S band group delays. The 
set B used X band group delays for an experiment 
with HOBART26 and X/S combinations for all other 
experiments. The set C used S band group delays. 
We got three series of 21 estimates of HOBART26 
positions. We formed the difference between esti- 
mates of HOBART26 vertical and horizontal coordi- 
nates between sets B and A, and between sets C 
and A. The differences for vertical component of 
site position estimates are given in Figure 6. 

The only difference between solutions in sets A, 
B, and C is the treatment of the ionosphere path 
delay. Therefore, we can consider the difference 
B — A and C — A are primarily due to the contribu- 
tion of unmodeled ionosphere path delay into esti- 
mates of site position. The average ratio of differ- 
ences {B — A)/{C—A) is 13.7 which is very close to 
the square of the ratio of effective ionosphere fre- 
quencies at X and S bands: (8.387/2. 250)^ « 13.9. 
This supports our argument that the differences 
C — A are due to the ionosphere path delay at 
S band. Since the differences B — A and C — A are 
nicely scaled with the square of effective frequency, 
we can extrapolate the estimates of the weighted 
root mean square (wrms) of site position time se- 
ries from solution C to L band by scaling them by 
(2.250/1.4)^ sa 2.6. If we were determining the po- 
sition of HOB ART 26 from One of these twenty- four 
experiment at L band, the wrms of position errors 
due to unaccounted ionosphere path delay would 
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Figure 6 The differences in the vertical coordinate esti- 
mate of station HOBART26 from a series of trial solutions, 
when position was estimated only from one experiment at 
a single baseline with PARKES. Hollow circles represent the 
difference between the X band only solution and the X/S 
dual band solution. The weighted root mean square (wrms) 
is 40 mm. Filled discs represent the difference between the 
S band only solution and the X/S dual band solution. The 
wrms is 490 mm. 

have been 1.3 ni for the vertical component and 
0.14 m for the horizontal component. We extrap- 
olate these numbers to our estimates of ASKAP-29 
and WARK12M positions. 

The differences in position estimates of WARK12M 
from VLBI and from GPS surveys, VLBI-GPS, 
are 1.66 m for the vertical component, 0.49 m for 
the east component, and —0.48 m for the north 
component. The difference in the vertical compo- 
nent can be explained by the effect of the iono- 
sphere. 

5 Summary and Future Work 

We obtained the first estimates of the positions of 
ASKAP-29 and WARK12M antenna reference points 
from VLBI observations. The random position la 
errors, 5-6 cm for the horizontal coordinates and 
20-40 cm for the vertical component, are close to 
that what one can expect from narrow-band ob- 
servations at L band. It is a pleasant surprise 
that the very first observations that followed the 
successful fringe test yielded a reasonable result. 
However, the systematic errors due to unaccounted 
ionosphere are significantly greater: 1.3 m for the 
vertical component and 0.14 m for the horizontal 
component. We also identified several problems 
with station equipment that will be fixed in the fu- 
ture. The use of a Rubidium frequency standard at 
ASKAP-29 caused a decorrelation at a level of 4-5% 
within 2 minute long scans. The frequency insta- 
bility at longer time intervals appeared negligible 
with respect to errors caused by the ionosphere. 

The existing L band (1.4 GHz) receiver on 
ASKAP-29 will be replaced in early 2011 with 
a new receiver comprising a purpose-built feed 



horn optimized for this antenna, with significantly 
higher efficiency and a noise calibration source. In 
2011 ASKAP-29 may be equipped with an X band 
(8.4 GHz) receiver, and the possibility of an up- 
grade to a higher frequency receiver is being con- 
sidered. 

In 2011, WARK12M was equipped with S/X 
dual-band receivers. Future geodetic experiments 
at these frequencies are planned for February 
2011. The participation of stations ASKAP-29 and 
WARK12M in VLBI experiments under absolute as- 
trometry and geodesy programs will permit an im- 
provement in the precision of results by two or- 
ders of magnitude and reach a millimeter level of 
accuracy. A more precise survey to provide 1-3 
mm accuracy of the tie vector between VLBI and 
GPS antenna reference points is planned. Com- 
bined with further VLBI observations it will then 
be possible to reconcile the currently observed dif- 
ferences between GPS and VLBI positions of the 
WARK12M antenna reference point. 

In November 2010, Telecom New Zealand made 
its 30-metre satellite Earth station available to 
AUT's Institute for Radio Astronomy and Space 
Research. It is a wheel-and-track beam-waveguide 
antenna built by NEC Corporation in 1984. It will 
be converted to a radio telescope capable of con- 
ducting both astronomical and geodetic research. 
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